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ABSTRACT 



Aims. To calculate transition rates from ground and excited states in neutral oxygen atoms due to electron collisions for non-UTE modelling of 
oxygen in late-type stellar atmospheres, thus enabling reliable interpretation of oxygen lines in stellar spectra. 

Methods. A 38-state /^-matrix calculation in Li' -coupling has been performed. Basis orbitals from the literature (Thomas et al.) are adopted, 
and a large set of configurations are included to obtain good representations of the target wavef unctions. Rate coefficients are calculated by 
averaging over a Maxwellian velocity distribution. 

Results. Estimates for the cross sections and rate coefficients are presented for transitions between the seven lowest LS states of neutral 
oxygen. The cross sections for excitation from the ground state compare well with existing experimental and recent theoretical results. 



Key words, atomic data 



1. Introduction 

Oxygen is one of the most abundant and important elements 
in the universe. Oxygen is widely accepted to be predomi- 
nantly produced in core-collapse supernovae resulting from 
massive stars, and thus detailed study of the chemical evolu- 
tion of oxygen allows us to understand star formation history 
and the initial mass function, as well as supernovae yields. 
The main observational evidence regarding the evolution of 
oxygen in the Galaxy comes from the measurement of oxy- 
gen abundances in the atmospheres of cool stars. The mag- 
nitude and trend of measured oxygen abundance with age or 
metallicity is not yet well constrained, in particular the re- 
sults at low metallicity differ between different studies (e.g. 
Israelian et al. 1998 Boesgaard et al.'l999 Nissen et al. l2U02l 
Garcia-Perez et al. 2006 1. The use of different spectral diagnos- 
tics is most likely the reason for these discrepancies. 

The use of the Oi triplet at 777 nm is appealling as the 
lines are strong and unblended, even in metal-poor stars. As 
they are in the visual region they are easily accessible with 
standard telescopes and spectrographs. On the other hand, it 
is well known that these lines are not formed in conditions 
of local thermodynamic equilibrium (LTE) in cool star atmo- 
spheres (see Kiselman 2001 and references therein). If the non- 
LTE formation of these lines could be understood, then they 
would be very useful diagnostics for measuring stellar oxygen 
abundances. However, we also note that the Unes originate from 
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highly excited levels, and are thus quite temperature sensitive, 
and furthermore may be saturated at higher abundances. 

The ability to solve the non-LTE problem is dependent on 
the availability of adequate atomic data for the radiative and 
collisional processes occuring to the oxygen atoms in the at- 
mosphere. The data for the collisional processes is the great- 
est stumbling block to reliable non-LTE modelling. Collisions 
with electrons are particularly important, due to their large fre- 
quency (due to high velocity) and efficiency (collisions are not 
adiabatic). Though collisions with hydrogen atoms are less fre- 
quent and expected to be less efficient (adiabatic) they may 
also be important due to sheer weight of number (typically 
A^H ^ lO'^A^,). Though Kiselman (2001 1 identifies the colU- 
sional cross-sections due to hydrogen as the most pressing issue 
since there is very little information on these processes (but see 
recent paper by Krems et al. 2006i, the situation for electrons 
is still far from ideal. Past studies of electron-impact excitation 
of neutral oxygen have predominantly focussed on excitation 
from the ground configuration. However, for the purposes of 
the non-LTE problem in stellar atmospheres, data for excited 
states are also needed. Previous work on the oxygen non-LTE 
problem have used either approximate general formulae (e.g. 
Allende Prieto et al. .2003) or a range of data calculated by dif- 
ferent methods (e.g. Carlsson & Judge 1993 Kiselman ll993> . 
In some cases data (except the simple formulae) do not exist, 
such as for the coupling between LS states of different spin, 
e.g. triplet-quintet system coupling. Kiselman performed test 
calculations for three stellar atmosphere models and found that 
intersystem collisional coupling was probably not important. 
Recent work by Fabbian et al. (private communication) con- 
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firms the results of Kiselman for the considered models, but 
finds that triplet-quintet intersystem coupling may be impor- 
tant in the case of metal-poor turnofF stars (e.g. Tf-g = 6500 K, 
log^ = 4, [Fe/H] = -3). 

In any case, it is clear that a consistently computed set of 
reliable data for electron-impact excitation of oxygen is highly 
desirable for non-LTE studies of oxygen line formation in stel- 
lar atmospheres. In this paper, we use the /^-matrix approach 
to calculate electron-impact excitation cross sections for low 
energy collisions for transitions between the seven lowest LS 
states of neutral oxygen, from which we calculate rate coeffi- 
cients for cool star atmosphere temperatures. We compare our 
results with experiment and those in literature, and discuss the 
differences. 



2. Calculations 

The electron scattering cross sections are calculated using the 
7?-matrix method (Burke et al. ll971t in LS -coupling with read- 
ily available computer codes. The internal region problem is 
solved using the the RMATRXl code (Berrington et al. ll995> . 
and the external region scattering problem is solved using the 
FARM code (Burke & Noble .1995,) . These papers, and refer- 
ences therein, should be consulted for detailed descriptions of 
the theory and the codes. We now describe the calculations in- 
cluding the choice of orbital functions, basis configurations, 
target states and calculation parameters. 

We employed the set of radial orbital functions from 
Thomas et al. (!1997> . This basis consists of Is, 2s and 2p or- 
bitals from Clementi & Roetti J1974> . and 3s, 4s, 3p, 4p, 3d 
and 4d spectroscopic orbitals calculated by Thomas et al.; this 
paper may be consulted for details. Their basis also included 
pseudo-orbitals labelled 5s, 5p and 5d; however, these orbitals 
are not employed in our final calculations, as will be discussed 
below. 

A well-documented problem (e.g. Plummer et al. 2004 
Berrington et al.T588 1 in /^-matrix calculations for moderately 
complex atoms such as oxygen, is the balance between a good 
description of the target wavefunctions, and the problem of 
pseudo-resonances arising due to an inconsistent description of 
the A^-electron target and (A^ -i- l)-electron problem. Computer 
memory limitations often lead to use of inconsistent sets of A^- 
electron and (A^ -i- l)-electron configurations. We note the exis- 
tence of methods which may allow this problem to be avoided 
or its effects drastically reduced, in particular the method devel- 
oped by Gorczyca et al. fl995l l and the B-spline /^-matrix ap- 
proach with non-orthogonal orbitals (e.g. Zatsarinny & Froese 
FischerUOOO Zatsarinny & Tayal 2001 1 a code for which was 
recently published by Zatsarinny (I2006t . With this in mind, 
we chose consistent A^- and (A^ -i- l)-electron configurations for 
our calculations. We included A^-electron configurations arising 
from single and double excitations from 2s^2p'*, but with max- 
imum occupations of orbitals for different states as given in 
Tabled For the (A^ + l)-electron configurations we account for 
all configurations arising from single, double and triple exci- 
tations from 2s^2p^ with the same maximum occupation num- 
bers as for the target states. 



Table 1. Limits on occupation numbers for orbitals. 



States 


3s 3p 


3d 4s 4p 4d 4f 


Max. Occupation Number 


2 2 


2 2 2 2 1 



As will be seen in the next section, our cross section re- 
sults are significantly different from those of Thomas et al., 
in particular we do not find the same slow rate of increase of 
the cross sections near the threshold. Thomas et al. chose con- 
figurations based on their expansion coeffients and the inde- 
pendent selection of particular configurations in the target and 
(A^ -I- l)-electron wavefunctions could have led to such pseudo 
resonances. To test this, we also performed calculations select- 
ing A^-electron configurations in the same manner as Thomas 
et al., i.e. based on their expansion coeffients, noting that the 
pseudo-orbitals 5s, 5p and 5d were included in this calculation. 
We obtained cross section results of similar form to those of 
Thomas et al., in particular we reproduced the slow rate of in- 
crease of the cross section near the threshold. An additional test 
calculation was performed in which the pseudo-orbitals were 
not used but the same configurations, except those involving 
the pseudo-orbitals, were selected. In this case the results were 
similar to our calculations where consistent A^- and (A^ -H 1)- 
electron configurations were employed. This suggests the in- 
clusion of the pseudo-orbitals may be the main cause of the 
slow increase near threshold; however, one must note that the 
degree of incompleteness between A^- and (A^-H l)-electron con- 
figurations is reduced due to the smaller set of basis orbitals. To 
resolve this definitively would require a calculation including 
the pseudo-orbitals and a large set of consistent A^- and (A^-H 1)- 
electron configurations, and this is presently beyond the com- 
puter resources available to us. Based on these considerations 
and results, we decided to limit our orbital basis to the spec- 
troscopic orbitals, enabling consistent sets of A^- and (A^ + 1)- 
electron configuration to be chosen. 

In the calculations we have included 38 target states in 
the close-coupling expansion which are listed in Table |2l The 
model includes 19 spectroscopic target states, and 19 additional 
eigenstates which were included to partially account for cou- 
pling to the continuum. Of these additional eigenstates, 14 cor- 
respond to observed autoionising states, while an additional 5 
(pseudo-) eigenstates are included so that at least one state of 
all symmetries which dipole couple to the lowest seven target 
states are included. We note that the inclusion of the autoion- 
ising and the pseudo-eigenstates does not affect the cross sec- 
tions greatly, and in particular the cross sections near threshold 
are basically unaffected. However, their inclusion lessens the 
magnitude of pseudo-resonances occuring above the ionisation 
threshold. Target energies for the 19 spectroscopic states are 
adjusted to the observed values in the scattering calculations. 
Calculations were also performed with reduced numbers of tar- 
get states. A calculation including only the 19 spectroscopic 
states, and a calculation including the 33 physical states (spec- 
troscopic and autoionising) were performed which we find use- 
ful for comparison below. 
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Table 2. The included target states and their calculated theoret- 
ical excitation energies and observed energies; the difference is 
also shown. The observed energies are taken from NIST (see 
text). 



State 


Energy 


Energy 


Energy difference 




(tlieory) 


(observed) (theory - observed) 




[Ryd] 


[Ryd] 


[Ryd] 




spectroscopic states 




2p*3p 


0.000 


0.000 


0.000 


Ip-* 'D 


0.154 


0.145 


-1-0.009 


2p'' 'S 


0.301 


0.308 


-0.007 


2p^^3s 


0.679 


0.672 


+0.007 


2p^^3s 'S" 


0.705 


0.700 


+0.005 


2p3 3p 5P 


0.795 


0.789 


+0.006 


2p3 3p ^"P 


0.820 


0.808 


+0.012 


2pMs ^S" 


0.875 


0.870 


+0.005 


2p34s ^S" 


0.882 


0.877 


+0.005 


2p33d ^D" 


0.893 


0.888 


+0.005 


2p33d 'D" 


0.893 


0.888 


+0.005 


2p34p 5p 


0.908 


0.903 


+0.005 


2pHp ''P 


0.919 


0.908 


+0.011 


2p33s ^D" 


0.933 


0.922 


+0.011 


2p33s 'D" 


0.947 


0.936 


+0.011 


2pHd ^D" 


0.942 


0.937 


+0.005 


2pHd 'D" 


0.943 


0.938 


+0.005 


2pHf 


0.943 


0.938 


+0.005 


2pHf 


0.943 


0.938 


+0.005 




autoionising states 




2p33s ^P" 


1.033 


1.039 


-0.006 


2p33s 'P" 


1.047 


1.057 


-0.010 


2p'4s ^D" 


1.135 


1.116 


+0.019 


2p'4s 'D° 


1.137 


1.120 


+0.017 


2p33d 'P" 


1.149 


1.125 


+0.024 


2p3 3d 'D" 


1.150 


1.133 


+0.017 


2p33d 'P" 


1.150 


1.133 


+0.017 


2p33d 'D" 


1.151 


1.133 


+0.018 


2p33d 'F" 


1.151 


1.134 


+0.017 


2p3 3d3S'' 


1.150 


1.134 


+0.016 


2s 2p^ 


1.197 


1.152 


+0.045 


2p34d 'P" 


1.199 


1.182 


+0.017 


2pHd 'F" 


1.200 


1.183 


+0.017 


2pHd ^P" 


1.220 


1.185 


+0.035 




pseudo-eigenstates 




3P 


1.086 






5P 


1.797 






5po 


1.898 








1.911 








1.930 







The R-matrix calculations are performed with 50 contin- 
uum orbitals for each of channel angular momentum of < Z < 
25. The iJ-matrix boundary is set at a = 56.4 aq. Partial waves 
with total angular momentum up to and including L = 20 are 
included, with full exchange included up to L = 12. Test cal- 
culations indicate that all cross sections are well converged at 
energies below 40 eV even with only L < 10 included. 

The accuracy of the wavefunctions may be tested, to some 
degree, by comparison with experimental observables such as 



excitation energies and oscillator strengths. In Table|2the the- 
oretical excitation energies are compared with experimental 
values taken from the NIST Atomic Spectra Database'. The 
results agree quite well with experiment, particularly for the 
lower lying states which are of interest here. In Table |3] the 
line strengths for transitions between low-lying states are com- 
pared with those from the NIST Atomic Spectra Database and 
Biemont et al. ( 1991 1. Again, the agreement is satisfactory, in- 
dicating that the wavefunctions are of acceptable quality. 

A coiTect treatment of long range dipole polarisation may 
be important for low energy electron-atom scattering calcula- 
tions. In the case of ions the long-range interaction is domi- 
nated by Coulomb forces, while for neutral atoms the lowest 
order interaction arises from the dipole moment induced in the 
atom by the electron's electric field. Accurate representation of 
the dipole polarisation requires that contributions from all cou- 
pled channels, including the continuum, be accounted for. In 
calculations, this may be achieved by the introduction of pseu- 
dostates, which mimic the contribution of a large number of 
channels (Damburg & Karule 1967 1. A correct representation 
of the dipole polarisability is very important for elastic scat- 
tering calculations, particularly the differential cross sections 
(e.g. Plummer et al. .2004. , Zatsarinny et al. 2006^. However, at 
least for oxygen, the available evidence points to this being less 
important in inelastic scattering processes. Tayal (2002i found 
that for the ^P-3s^S° and ^f-3s^D° transitions, inclusion of 
coupling to the continuum reduces the cross sections by 5 to 
15% and 5 to 27% respectively. Zatsarinny & Tayal (2001) and 
Plummer et al. (2004) found small effects (< 10%) for the ^P- 
'D, and ^S-^D transitions. Tayal J2004> concludes that 

the effect of coupling to the continuum is greater on resonance 
transitions than on forbidden transitions. The effect is also ex- 
pected to be mostly important above the ionisation threshold 
(e.g. Tayal 200? I, the differences seeming to be largest near the 
maximum in the cross section, but negligible within a few eV 
of the threshold. 

In our calculations we have not included any polarisation 
pseudostates, and only a fraction of the true dipole polaris- 
ability for each state is included. For example, for the ground 
state the experimental value for the static dipole polarisability 
is found to be a = 5.4 + 0.4 a.u. (Alpher & White ll959> . while 
theoretical values are typically between 4.6-5.4 a.u. (see results 
collected in Thomas et al. I1997> . We calculate for our model 
a - 1 .68, and thus only about a third of the polarisability of this 
state is included. Based on the above, noting that the majority 
of the transitions are forbidden, we expect that this may lead to 
an overestimate of cross sections for allowed transitions by of 
order 10%. Given that we are interested in quite low tempera- 
tures and thus low collision energies, the cross sections within 
a few eV of the threshold will dominate in the rate coefficients 
and thus the errors would most likely be even smaller. 



' National Institute of Standards and Technology, Atomic Spectra 
Database,_http://physics.nist.gov_ 
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Table 3. Comparison of computed line strengths with data from 
the NIST compilation, and where available the theoretical cal- 
culation of Biemont et al. (I1991> . 





Transition 


5 

(this work) 


5 

(NIST) 


S 

(Biemont et al.) 




-2p^3s^S° 


2.211 


2.001 




2p33s ^S" 


— 2p33p^P 


130.8 


128.4 


123.7 


2p33s ^S" 


— 2p33p3p 


81.8 


86.3 


89.5 



3. Results and Discussion 

The computed cross sections are plotted in Fig. ^ Note, that 
in both Figs, [l] and |2l some regions which were affected by 
pseudo-resonances at energies above the ionisation threshold 
were either smoothed or interpolated to remove any sharp tran- 
sient behaviour which is unlikely to be physical. The results of 
the 38-state and 19-state calculations are shown, both smoothed 
in the same manner. Comparison of these results, along with 
the 33-state calculation, allows pseudo-resonances to be iden- 
tified as features which change between the calculations thus 
are very likely to be spurious (due to incompleteness of the 
model), while features that remain unchanged are most likely 
real. In general, we may note that the near threshold behaviour, 
most important for the temperatures of interest, are not greatly 
affected. However, it is clear that many of the cross sections 
are probably affected by pseudo-resonances particularly at to- 
tal energies around 15-20 eV. Note, the feature at around 14 eV 
in the ^Z)-'S cross section is a real feature, it is a resonance in 
the partial cross section (see e.g. Thomas et al. ll997> . 

As expected, the singlet-quintet transitions have null 
cross-sections since in LS' -coupling total electronic spin 
must be conserved. In addition, the ^S—^^S" transition is 
forbidden as transitions are forbidden in general 

(Goddard et aL lTWlT l. 

Experimental results at low energy are available for some 
transitions from the ground state. We compare our results with 
experiment in Fig.|2l In the case of the ^P—>^S'' transition, we 
also plot the cross section from the 4 1 -state 7?-matrix with pseu- 
dostates calculation tabulated in Johnson et al ( 2003aJ and de- 
scribed in Tayal 12002 1. 

In general the computed data compare quite well with the 
experimental data considering the experimental error bars. For 
the — >'D transition, the experimental results suggest a peak 
in the cross section at an impact energy of about 6 eV. A sim- 
ilar peak is seen in the computed cross section, though only 
roughly half the magnitude. Compared to the cross section of 
Thomas et al., the peak is about of the same magnitude but 
at lower energy, with our cross section not showing the slow 
rate of increase in the cross section of that of Thomas et al., 
but rather a sharp increase at the threshold. The cases of the 
^f— >'5 and 'D— (see Fig. [3 transitions are quite similar, 
again rather than the slow rate of increase in the cross section 
just above threshold of Thomas et al., we find a sharp increase. 

The cross section for the ^P—^^S" transition compares well 
with experiment at low energy, but is significantly larger than 



the experimental values at energies at and above 20 eV. The 
cross sections also compare well with those of Zatsarinny 
& Tayal (2002"), having the same magnitude and form near 
threshold. The peak at about 18 eV seems to be the result 
of pseudo-resonances, as judged from the fact that the peak 
varied significantly between the 19-, 33- and 38-state calcula- 
tions (see Fig. We note the peak is not present in the re- 
sults of Zatsarinny & Tayal, which uses the non-orthogonal or- 
bital method which should be more robust against such pseudo- 
resonances. Similarly, the ^P—>^P and ^P—>3p^P results seem 
to be affected by pseudo-resonances at total energies of about 
18 eV. However, near threhold the results are similar to those 
of Zatsarinny & Tayal. 

The resonance transition ^P—>^S ° is the best studied, with a 
significant number of experimental and theoretical studies. As 
seen in Fig. |2] our result compares well with the experimen- 
tal results, both in terms of magnitude and form of the cross 
sections. The cross section compares well with the theoretical 
result of Tayal (2002i at energies between 15 and 30 eV, which 
is also plotted. Our cross sections are around 5-10% larger, as 
expected (see discussion in previous section) due to our neglect 
of a large part of the coupling to the continuum. Regarding the 
near threshold behaviour, our results are very similar to those 
of Zatsarinny & Tayal (2002i, the cross sections having very 
similar magnitude and form, including a sharp resonance at the 
threshold. 

Rate coefficients are computed by integrating a Maxwellian 
velocity distribution with the cross-sections, and are presented 
in Table |4] The calculations are done assuming LS coupling, 
and as such no fine structure is included. Fine structure of 
certain levels of interest are often included in non-LTE mod- 
elling, and one may wish to somehow redistribute the rate co- 
efficent between fine structure levels. It is often the practice 
to set fine structure transition rates to large values such that 
the fine structure levels have populations according to their sta- 
tistical weights (e.g. Kiselman 1993). As such, it is irrelevant 
how the rates are distributed so long as the total rate is con- 
served. However, we note that even neglecting relativistic ef- 
fects in the Hamiltonian, a simple redistribution of the rate co- 
efficients is not possible as a transformation from LS to Jaj 
coupling (where and j are the initial total angular momenta 
of the atom and scattered electron respectively) must be done 
for each partial wave and scattered electron angular momentum 
(e.g. Saraph frgTot . 

It is of interest to compare these results, at least in broad 
terms, with those that have been used in past non-LTE calcu- 
lations. Below we make comparisons at T - 5000 K, which 
is representative of the typical temperatures of interest in cool 
stellar atmospheres. First we make a comparison with the pre- 
dictions of approximate yet general and simple formulae used 
by Allende Prieto et al. (2003 1; this paper may be consulted 
for the formulae. For allowed radiative transitions they used a 
formula due to Mihalas (1972i, which has been derived from 
the approximate formula of van Regemorter J1962> adopting 
an effective Gaunt coefficient of 0.25. For forbidden transi- 
tions an effective collision strength of T,y = 0.05 (for exci- 
tation) has been arbitrarily adopted. The collision rate coeffi- 
cient is computed from the effective collision strength via what 
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Fig. 1. The cross sections for electron-impact excitation as a function of initial electron energy between the seven lowest target 
states. The full lines are the 38-state calculation, while the dotted lines are the 19-state calculation. The cross sections are in units 
of 10-1^ cml 



Allende Prieto et al. have referred to as the Eissner-Seaton for- 
mula (Eissner & Seaton n"974> . We calculated rate coefficients 
at 5000 K, and these along with the ratio of these values to our 
values are presented in Table|5l We see for allowed transitions, 
the ratios range from about 0.6 to 3.0. As mentioned, in LS- 
coupling transitions between singlet and quintet states, and the 
'5 — >^ S° transition, are not collisionally allowed and thus we 
cannot make a comparison here. For the remaining forbidden 
radiative transitions, the ratios range between roughly 0.003 
and 60 indicating much poorer agreement in general. Based 
on the evidence here, it appears that the approximate formula 
based on van Regemorter's formula gives at best order of mag- 
nitude estimates, while the simple approximation based on an 
arbitrarily chosen constant T,j, is considerably less reliable. 

We also compared with the data in the model atom prepared 
by Carlsson & Judge ( 1993 1 at similar temperatures. The first 
notable aspect of the comparison is that for several transitions 
no collisional coupling was included due to lack of data. For 
the transitions ^P—>^D, '5 and 'D— they used data from a 
compilation by Mendoza lil981) , which in fact originate from 
Berrington & Burke (I1981> . At 5000 K, the corresponding ex- 



citation rate coefficients are 1.82 x 10"", 1.27 x 10"'"^ and 
1.04 X 10"" cm s"' respectively. The ratios with respect to 
our values are roughly 1.3, 0.4 and 0.8 respectively. For the 
^P—>^S°, transitions, they calculated collision rates from 
the cross sections of Rountree J1977> . The cross sections of 
Rountree are significantly larger at the threshold than modern 
calculations, including our results, resulting in their rates being 
a factor of about 2.5 greater than ours. For the ^5°— and 
^S"—^^P transitions, Carlsson & Judge computed rates using 
the impact parameter method of Seaton (I1962> . The rates are 
about a factor of five and two greater than ours, respectively. 

Zatsarinny & Tayal ('2003"l have published effective col- 
lision strengths based on calculations using the B-spline R- 
matrix approach with non-orthoganal orbitals. They provide 
data for transitions between the ^P, 'D and '5 states of the 
ground configuration and from these states to 23 other excited 
states. We compare our results with theirs at 5000 K in Table|6l 
The data are in good agreement, agreeing within a factor of 3 
at worst, often much better, the mean ratio (o'v)zT/(o'i')thiswork 
being 1.23 with a standard deviation of 0.67. 
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Fig. 2. The collision cross sections as a function of initial electron energy compared with experiment. The full lines are the 
results of the present work. Experimental and other theoretical results in each case are from: (a) Doering ( 1992), circles; Shyn 
& Sharp ( 1986 squares; in the latter case, the error bars are slightly offset for clarity, (b) Doering & Gulcicek ( 1989ai, circles; 
Shyn et al. ( 1986 1, squares; in the latter case the error bars are slightly offset for clarity, (c) Doering & Gulcicek ( 1989b i, circles, 
(d) Doering & Yang (2001 1, circles; Johnson et al ( 2003a i, small filled squares; model-based results from experimental electron- 
impact-induced emission cross sections from Johnson et al (pOOSbi, open squares; 41-state calculation of Tayal J2002> . thick 
dotted Une. (e) Gulcicek et al. fWgSt . circles, (f) Gulcicek et al. al988l . circles. 



Estimating the errors in the computed rate coefficients is 
difficult. The generally good agreement of our cross section 
results with experiment, particularly at low energy, would indi- 
cate that our rate coefficients are of quite acceptable accuracy. 
Errors due to neglect of coupling to the continuum and effects 



of pseudo-resonances at energies above the ionisation thresh- 
old are expected to be small as these energies hardly contribute 
to the rate coefficient for the temperatures of interest. Perhaps 
the best estimate of the error is obtained from the scatter among 
different /?-matrix calculations where available. From the com- 
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Table 4. Rate coefficients (crv), in units of cm s ' , for electron-impact excitation of neutral oxygen, for selected temperatures in 
the range T = 1000-50000 K. Note, a{b) = a x 10^ 



Initial 






Final state 






state 


2p'' 'O 


2p4 'S 


38^5° 3s ^5° 


3p V 


3p3/' 



1000 K 



2p'*3/' 


5.72(- 


-21) 


5.13(- 


31) 


1.03(- 


55) 3.40(- 


-57) 


6.35(- 


-64) 


5.48( 


-65) 


2p'* '£» 






1.91(- 


20) 




0. 3.34(- 


-49) 




0. 


6.34( 


-56) 


2p'* '5 












0. 


0. 




0. 


8.85( 


-45) 


38^5° 












9.98(- 


-10) 


I.IK- 


-15) 


1.98( 


-17) 


3335° 
















9.91(- 


-14) 


3.44( 


-15) 


3p^P 












3000 K 








2.56{ 


-09) 


2p'* V 


2.98(- 


-13) 


4.28(- 


17) 


9.33(- 


25) 3.46(- 


-25) 


9.03(- 


-28) 


4.89( 


-28) 


2p'* '£) 






4.17(- 


13) 




0. 6.87(- 


-24) 




0. 


9.66( 


-26) 


2p'* '5 












0. 


0. 




0. 


4.53( 


-22) 


38^5° 












2.36(- 


-08) 


6.23(- 


-10) 


4.43( 


-11) 


3s ^S" 
















1.65(- 


-09) 


9.26( 


-10) 


3p5p 












5000 K 








3.51( 


-08) 


2p'*3/' 


1.38(- 


-11) 


3.07(- 


14) 


1.30(- 


18) 7.53(- 


-19) 


1.53(- 


-20) 


1.26( 


-20) 


2p'' '£» 






1.24(- 


11) 




0. 6.11(- 


-19) 




0. 


9.93( 


-20) 


2p^'5 












0. 


0. 




0. 


1.41( 


-17) 


38^5° 












3.89(- 


-08) 


1.02(- 


-08) 


8.34( 


-10) 


38 35° 
















1.09(- 


-08) 


1.29( 


-08) 


3p^/' 












8000 K 








5.83( 


-08) 


2p'* V 


1.24(- 


-10) 


1.34(- 


12) 


3.42(- 


15) 2.66(- 


-15) 


1.72(- 


-16) 


1.88( 


-16) 


2p'* '£» 






8.25(- 


11) 




0. 3.75(- 


-16) 




0. 


2.63( 


-16) 


2p'* '5 












0. 


0. 




0. 


7.98( 


-15) 


38^5" 












4.57(- 


-08) 


4.97(- 


-08) 


3.95( 


-09) 


38 35° 
















2.86(- 


-08) 


5.65( 


-08) 


3p'P 












12000 K 








6.98( 


-08) 


2p'*3/> 


4.31(- 


-10) 


1.10(- 


11) 


2.53(- 


13) 2.47(- 


-13) 


2.99(- 


-14) 


3.96( 


-14) 


2p'* '£» 






2.31(- 


10) 




0. 1.50(- 


-14) 




0. 


2.18( 


-14) 


2p'* '5 












0. 


0. 




0. 


3.30( 


-13) 


38^5° 












4.48(- 


-08) 


1.20(- 


-07) 


8.36( 


-09) 


38 35° 
















4.38(- 


-08) 


1.26( 


-07) 


3p5p 












20000 K 








6.89{ 


-08) 


2p'*3/> 


1.18(- 


-09) 


5.84(- 


11) 


7.37(- 


12) 9.56(- 


-12) 


1.91(- 


-12) 


3.00( 


-12) 


2p'* '£) 






5.15(- 


10) 




0. 3.24(- 


-13) 




0. 


8.02( 


-13) 


2p'* '5 












0. 


0. 




0. 


6.05( 


-12) 


38^5° 












3.74(- 


-08) 


2.38(- 


-07) 


1.30( 


-08) 


38^5° 
















5.25(- 


-08) 


2.31( 


-07) 


3p^/' 












50000 K 








5.81{ 


-08) 


2p'*3/> 


2.71(- 


-09) 


2.35(- 


10) 


1.32(- 


10) 3.04(- 


-10) 


8.95(- 


-11) 


1.77( 


-10) 


2p'* '£) 






1.05(- 


09) 




0. 5.85(- 


-12) 




0. 


2.57( 


-11) 


2p'* '5 












0. 


0. 




0. 


6.50( 


-11) 


38^5° 












2.03(- 


-08) 


4.12(- 


-07) 


1.26( 


-08) 


38^5° 
















4.02(- 


-08) 


3.62( 


-07) 


3p^P 




















3.23( 


-08) 



parison with the rates of Zatsarinny & Tayal J2003t and the 
rates of Berrington & Burke ('1981'i we would estimate that the 
rate coefficients have errors of the order of 70% at the temper- 
atures of interest. 



It is difficult to predict the impact that these new data 
will have on the results of non-LTE calculations in stellar at- 
mospheres. The impact on abundance analysis of solar-type 
stars, including metal-poor stars, is presently being investi- 
gated by Fabbian et al. (in preparation). It would also be of 
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Table 5. Rate coefficients (crv) in units of cm^ s ' for electron-impact excitation of neutral oxygen for T - 5000 K calculated 
using approximate formulae following Allende Prieto et al. J2003> . and the ratios of these rate coefficients to our rate coefficients 
in TableEl«crv)3pp,.ox/(o-v)this„ork)- Note, a{b) = a x 10*. 



Initial 






Final state 






state 




'0 


2p^ 'S Ss^S" 3s 


3p V 


3p V 








(crv) from approximate formulae 




2p* V 


7.04(- 


12) 


4.04(-14) 4.06(-19) 4.91(-19) 


1.00(-20) 


5.62(-21) 


2p'* '£) 






7.01(-12) 7.04(-17) 2.94(-17) 


1.74(-18) 


9.75(-19) 


2p'* '5 






6.13(-14) 2.56(-14) 


1.51(-15) 


8.50(-16) 


38^5° 






5.11(-10) 


2.18(-08) 


1.69(-11) 


3s 








1.20(-10) 


3.79(-08) 


3p V 










2.29(-10) 








{o'v)approx /(fv)this work 








5.10(- 


01) 


1.28(+00) 3.10(-01) 6.48(-01) 


6.48(-01) 


4.43(-01) 


2p'* '£» 






5.53(-01) - 4.79(+01) 




9.77(-i-00) 


2p'* '5 










5.98(-i-01) 


3s ^5" 






1.31(-02) 


2.14(+00) 


2.02(-02) 










1.09(-02) 


2.94(-i-00) 



3p 3.92(-03) 



Table 6. Rate coefficients (crv) in units of cm-' s ' for electron-impact excitation of neutral oxygen for T = 5000 K from 
Zatsarinny & Tayal ( 2003i) for transitions from the ground configuration, and the ratios of these rate coefficients to our rate 
coefficients in Table|2({crv)zT/(crv)thiswork)- Note, a(b) - ax 10*. 



Initial 
state 


2p'' '0 


2p'* '5 


Final state 


2,p^P 


3p V 








{crv) from Zatsarinny & Tayal 




2p'*V 


2.05(-ll) 


1.16(-14) 


1.94(-18) 1.15(-18) 


3.42(-20) 


2.86(-20) 


2p'* '£) 




7.79(-12) 


0. 3.45(-19) 


0. 


7.29(-20) 


2p'* '5 






0. 0. 


0. 


1.63(-17) 








<0"V>ZT/<0"V>thiswork 






2p'*3/' 


1.49(-h00) 


3.67(-01) 


1.48(+00) 1.52(+00) 


2.21(+00) 


2.25(+00) 


2p'* '£» 




6.15(-01) 


- 5.62(-01) 




7.31(-01) 


2p'* '5 










1.14(-i-00) 



interest to investigate the impact on predicted solar centre-to- 
limb behaviour of the oxygen triplet lines (see e.g. Allende 
Prieto et al. 2004). 

We have presented estimates of cross sections and rate coef- 
ficients for all transitions between the seven lowest LS states of 
neutral oxygen. It should be noted that the calculations produce 
data for transitions between all target states listed in Table |3 
We restricted the data presented here to the lowest seven states 
for two reasons. First, the quality of the results for transitions 
involving higher states is certainly lower. In particular, the ef- 
fect of pseudo-resonances on the rate coefficients would be ex- 
pected to be larger, since the thresholds are nearer the ionisa- 
tion energy. Secondly, the transition rates among the low-lying 
states are expected to be most important for the non-LTE prob- 
lem of interest, in particular the modelling of the 777 nm triplet 
3s ^5" ^ 3p ^P. However, the computed rates for the transi- 
tions involving the higher states should give at least order of 
magnitude estimates and may be obtained on request from the 
author. 



More calculations and experimental data for electron- 
impact excitation of oxygen would be valuable, particularly 
from excited states where few or no other comparisons are 
available. In particular calculations for transitions between ex- 
cited states using the Z?-spline 7?-matrix method with non- 
orthogonal orbitals, with more careful consideration of short 
range correlation and long range polarisation effects would be 
important. Calculations for excitation by hydrogen collisions 
are also needed. Recent work by Krems et al. ( 2006 ) calculated 
rates for a single transition, but again as for electrons, for stellar 
atmospheres applications, data are needed for more low-lying 
states. 
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